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Abstract 

The entropy production and the variational functional of a Laplacian diffusional field around the first 
four fractal iterations of a linear self-similar tree (von Koch curve) is studied analytically and detailed 
predictions are stated. In a next stage, these predictions are confronted with results from numerical 
resolution of the Laplace equation by means of Finite Elements computations. After a brief review of 
the existing results, the range of distances near the geometric irregularity, the so-called ’’Near Field”, a 
situation never studied in the past, is treated exhaustively. We notice here that in the Near Field, the usual 
notion of the active zone approximation introduced by Sapoval et al.^’^ is strictly inapplicable. The basic 
new result is that the validity of the active-zone approximation based on irreversible thermodynamics is 
confirmed in this limit, and this implies a new interpretation of this notion for Laplacian diffusional fields. 
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I. INTRODUCTION 


On the past decade, an important field of interdisciplinary research has been developed based 
on new geometry concepts, namely ’’the fractal” geometry of nature, pioneered by Benoit B. 
Mandelbrot. The fractal geometry approximates well enough the naturally discovered disordered 
morphologies. For instance, the terminal part of the respiratory system of mammals, biological 
membranes, porous electrodes or catalysts, provide some characteristic examples^”^^. 

This implies a direct link to Laplacian transfer, as one could mention that under ordinary 
conditions many familiar transport phenomena such as diffusion and heat conduction are described 
in the steady state by Laplacian fields. The Laplace equation is a linear equation, however the 
complexity of the problem is due to the geometrical irregularities of the boundaries. Due to this 
fact, an important number of studies on the distribution of Laplacian fields around geometrically 
irregular - and eventually fractal - objects, have been justified^’^^“^^. 

An old mathematical conjecture states that the Laplacian field around a fractal object in two 
dimensions is itself of multifractal nature. This conjecture, after Makarov’s theorem is a mathe¬ 
matical fact. In a paper considered classical by now, B.B. Mandelbrot explored for the first time 
in detail the Laplacian field around a deterministic self-similar tree with his famous ’’zebroide” 
(logarithmic) representation^^. This work was the starting point for many authors who examined 
not only diffusion but also reaction-diffusion systems^’^^“^^. 

Also, specialists from mathematics (harmonic analysis) presented landmark studies on the dis¬ 
tribution of harmonic measure around irregular boundaries. In particular, N.G. Makarov proved 
a theorem stating that, whatever the shape of an irregular (simply connected) boundary in two 
dimensions might be, the active zone in which most of the flux generated by a Laplacian field 
is concentrated, scales as a length^^“^^. More importantly, this result holds exactly also for the 
’’mathematical” fractal object. In the ’’multifractal” picture of the Laplacian field this amounts to 
the equality of the first multifractal exponent (the ’’information dimension”) exactly to the value 
1. Later on, B. Sapoval proposed another interpretation of Makarov’s theorem on the grounds of 
’’information ensemble”Moreover, the active zone approximation tells us simply that some 
parts of an irregular surface are not really contributing, remaining ’’passive”, so that all thermo¬ 
dynamic properties can be estimated by considering only the ’’active zone”. The notion of the 
active zone received ample experimental support through electrochemistry experiments leading to 
its visualization with electrodeposition of copper^^ (optical absorption). According to the theorems 
mentioned above, for a Laplacian field the active zone is always a good approximation. J. Bourgain 
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generalized this result in higher dimensions^, but with a scaling exponent which depends on the 
dimensionality of space and is not necessarily an integer. 

Continuing in the same lines, Makarov and Jones proved that the harmonic measure around 
some self-similar objects is multifractal. Recently, due mainly to work by Sapoval’s group, the 
problem of the multifractal spectrum around von Koch’s curve in two (and then in three) dimensions 
has been reconsidered on the grounds of hierarchical random walkers. Important evidence that 
in three dimensions there is a nontrivial behavior of the information measure and a non-integer 
dimensionality equal to 2.007 is given by Grebenkov et al.^^. 

In32 

an alternative way to characterize the complexity of Laplacian transport across irregular 
boundaries was proposed, based on irreversible thermodynamics. More specifically, the focus was 
set on the dissipation generated by the underlying process, as the irregularity of the boundary 
is increased. Dissipation is measured by the entropy production, arguably the central quantity 
of irreversible thermodynamics. The correct measure of complexity in this case is the entropy 
production per unit surface. In that study, as well as in^^, attention was followed to the ’’Far 
Field”, that is in distances from the object which far exceed the length of the minimal irregularity 
of the surface (the cut-off threshold). 

The main purpose of this work is firstly to revisit and generalize the results of^^ and secondly 
to explore the behavior of the entropy production functional in the region of the ’’Near Field”, i.e. 
for distances between the membranes (or electrodes) which are smaller than the smallest detail 
of the prefractal iteration. It should also be noticed that the chart of the spatial distribution of 
the local entropy production for an increasingly irregular object in the Near Field does not exist 
anywhere in the literature to the best of our knowledge. 

The paper is organized as follows. General considerations and detailed predictions for the active 
zone concept of an entirely diffusional field are presented in Sec. II. In Sec. Ill the near domain 
predictions are examined in detail and in Sec. IV these predictions are corroborated by numerical 
results obtained with Finite Elements and the chart of the local entropy production is portrayed. 
Finally, in Sec. V we draw the main conclusions and we discuss some more general issues. 

II. THE ACTIVE ZONE OF DIFFUSION FIELDS: GENERAL CONSIDERATIONS FOR 
THE ’’FAR” AND ’’NEAR DOMAIN” 

The simplest setting in which possible effects of complex geometry on entropy production can 
be identified is given in Figure 1. A concentration difference ci — C 2 is applied across the vertical 
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boundaries of a cell, over a characteristic length Ly. The cell obeys to zero flux boundary conditions 
along the horizontal direction x, but there is now a geometric irregularity consisting of extending 
the horizontal characteristic length Lx by a bump in the middle. 

Under the condition that c satisfies in the steady state Laplace equation 

V2c = 0, (1) 


it was shown^^ that, on the grounds of nonequilibrium thermodynamics, the full expression for the 
entropy production reads 


P^Dk 


{Vcf 


dr. 


y 

while the variational functional which is extremal in the steady state, is given by 


( 2 ) 


Q = J (Vc)^*. (3) 

y 

Before addressing the role of complex boundaries in the entropy production and the variational 
functional, we evaluate both c and P in the simple reference case of a two-dimensional box of 
height Ly and length Lx- The corresponding solution reads 


. . Cl — C2 .. X 

c(y) = —y + C 2 , (4) 

Ijy 

where ci = ci{Ly) and C 2 = C2(0). In addition, the entropy production can be calculated exactly, 
yielding 


P = - Cl) In —, 

Ly C\ 

whereas the variational functional in this case reads 


(5) 


Q = L(c2-ci)^. (6) 

Ijy 

For the setting of Figure 1, we have shown that in the ’’Far field” the entropy production would 
keep the structure of eq.(5) as far as y-dependence goes, but the proportionality factor Lx would 
be modulated based on Makarov’s theorem by an additional size-independent factor A depending 
on geometry and on concentrations only 


I C2 

F = DkA(ci,C 2 ;g)L^ — (c 2 - ci) In—. 

Ly Cl 

Under the same conditions the variational functional Q takes the form 


(7) 


Q = B{g)Lx — {c2 - cif, 

-L/^i 


( 8 ) 
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where now the constant B depends on the geometry only. It should be emphasized however, that 
the arguments supporting the two equations above are mainly heuristic and based on the grounds 
of irreversible thermodynamics and dimensional analysis. It is perhaps interesting to notice that 
a formal proof of these equations is lacking both in terms of random walks on a lattice or in 
continuous time random walks. 

In addition to the above toy-model fractal geometry, one might study more complex objects, 
i.e. higher fractal generations. As an illustration we present in Figure 2 the second and the fourth 
fractal iterations of the fractal generator of Figure 1. These irregular objects will play the role of 
counter-membranes in the sequel. 


C 
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FIG. 1: Schematic representation of the diffusion cell for the first fractal iteration of a fractal generator. 
The dimensions of the cell, = 3/ and Ly = 21, as well as the boundary conditions are indicated. For this 
fractal Dp = 

We now focus on the case of mild diffusion. The purpose is to estimate the coefficients A and 
B analytically in this limit. To this end we follow the so-called ’’independent field approximation” 
(IFA). This is a coarse-graining argument based on compartmentalization of the full continuous 
space in which diffusion takes place into a finite number of properly selected non-overlapping and 
non-interacting rectangular regions. Here, one can write down closed expressions for the entropy 
production and the variational functional assuming linear concentration profiles. The nonlinear 
dependence of the field in terms of the spatial coordinates arising from the irregularities of the 
boundaries is thus approximated by piecewise linear functions, entailing discontinuities of the 
equipotential lines at the boundaries between the cells. 
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FIG. 2: Shown are (a) the second and (b) the fourth fractal iterations of the fractal generator of Figure 1. 



FIG. 3: Schematic representation of the diffusion cell for the first fractal iteration of a fractal generator. The 
distance Ly = ml (here / = 1) is successively reduced from left to right: a) Ly = m = 0.3 to b) Ly = m = 0.1 
and c) Ly = m ^ 0. 


Let us illustrate how the IFA works for the first few fractal generations of the cell. We separate 
the cell into the three rectangles shown in Figure 3. Accepting linear concentration profiles in 
each of the side parts (ai) and (as) and applying eq.(5) with Lx = Ly — ml (i.e. setting 
F = Dk{c 2 — Cl) In one finds for the respective entropy production 


Pai — Pas — 

m 


( 9 ) 


Assuming furthermore a linear half penetration inside the pore (that is, a linear concentration 
profile until the middle of the pore^^“^^, thereby neglecting the remaining active zone) the central 
region (a 2 ) from eq.(5) with Lx — I and Ly — ml + ^ reads 


ml + ^ + 1 

Thus, the total entropy production of the cell becomes 


( 10 ) 


1 

Ptot — Pai + Pa2 + Pas — 2F" — . x • (11) 

m(2m + 1) 

Applying now eqs.(5) and (6) for the whole cell with Lx — 3/ and Ly — ml we find (remember that 
we deal here with mild diffusion, see eq.(7)) that 


m 


(12) 
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So that, for the first iteration (Figure 1) the coefficient B is of the form 

6m + 2 2 2 m 

3(2m + 1) 3 3 2m + 1 

Similarly, for the second fractal iteration (Figure 2a) we obtain 

324m2 + 150m + 12 4 2 T 1 

= -TTT-^ ^ 7: + ^rn\ - -- + 


and for the third 


9(6m + l)(6m + 3) 9 3 [2m+ 1 1 + 6m 




1 2 4 ■ 

2m + 1 6m + 1 18m + 1 ’ 


where the superscript (i) denotes the i — th fractal iteration. Proceeding in the same lines one can 
find and prove by induction that for the i — th fractal iteration it holds 

k=l 

The superscript (i) denotes the corresponding fractal generation (i.e. i = 1 is for the first fractal 

generation, i = 2 is for the second etc) and m indicates the respective inter-electrode distance. 

(i) 

Pausing for a moment the analysis of •> we investigate some features for the total entropy 
production of the fractal generations. Especially, one can easily observe that for the first few fractal 
iterations the total entropy production is 

Pill = Pill = ^Pif + ‘^Pif + Pif^ ( 17 ) 


Pill = 8Pif) + 4Pif + 2P^f + Piy Pltl = 16Pi^) + 8Pi^) + AP^f + 2Pif + pP\ (18) 

where the superscript {%) refers to the fractal iteration and the index ai corresponds to the geomet¬ 
rical structure that one can encounter in each fractal (see also Figure 3). For example, the index 
ai (or equivalently aa) refers to the left (right) side part of the first fractal generation, while a 2 
denotes the central region. Similarly, the indices a 4 , refer to the central regions of the second 
and third fractal generations respectively etc. The above procedure implies that for the i — th 
fractal iteration the total entropy production should be written as 

pS = 2*P« + 2*-ip« + 2*-2p« + 2*-3p« + (19) 

or in a more compact form 


E 

k=0li—k^0 


2*“^ py. 
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Here, the superscript (i) denotes the corresponding fractal generation (i.e. i = 1 for the first fractal 
generation, i = 2 for the second etc), whereas k corresponds to each geometrical structure that 
appears in the corresponding fractal generation as explained above. Moreover, we can find that 
the entropy production for each geometrical structure can be written as 


p(^ = - J— and Pjf'^ =-- 

1 32 k 2-3^ +1 


, A; ^ 1, i ^ A;, 


( 21 ) 


while the above polynomials are connected through the recurrence relations 




_ -I I 4-3^m J > k k > ^ and 1 > h h > ^ (22) 

p(l) - 3^+i+2-3^m’ ^ ^ ^ ^ - 3fe+2-3^m ^ ^ ^ 1 . 


Note also that the total entropy production for each fractal generation and geometrical structure 
can be written as 




-)k—l 


k=l 


+ 2 • 


(23) 


(i) 

Coming back to the investigation of the coefficients we observe that the first derivative with 
respect to m is positive 


dm 


-)k—l 


iy- _ 

3 ^(2.31,-1+ i) 


m 


2 ’ 


(24) 


implying that is an increasing function of m. We should also observe here that for m= 1, 2, 
3 one finds again the results of^^, (geometries for the Far field) 


- 0 . 8889 , Bf^ = ^ ~ 0 . 8571 ,...,= Itatll'llt - 0 . 8471 , 
^9 ^7 ’ ’ 1 289 , 739,835 

with a limiting value (in the context of the IFA approximation) 

(10) ^ 22, 614,533,817,495, 831, 212,892,038 ^ 

^ 26, 698,620, 787,529, 312,076,958,125 


(25) 


(26) 


III. PREDICTIONS FOR THE NEAR FIELD 


As also pointed out in the introduction, the ’’Near Field” is a situation not studied before in the 
literature. Supposing that in the phenomenon there is no new physics, we can apply our real space 
renormalization scheme (IFA) in order to obtain predictions and compare them with the finite 
element predictions, considered here as reference values. More specifically, we can easily confirm 
that for m — 0.3 the predictions of our real space renormalization scheme are 


^ 0.7916, = 7^ - 0.7123, ...,B^^l = 7^77^ - 0.6821, 

0.3 24 504 9069732 


(27) 



approaching asymptotically the value (in the framework of the IFA) 

„(io) _ 2,574,446,468,051,339,779,667,687 ^ _ 

3,776,514,354,889,019, 285,033,772 


Similarly, for m=0.1 the predictions of the renormalization are 


r(1) - 

^0.1 - ig 


0 . 7222 , 5 ^ = ^ 


('5') 37388 

0.5833,.= —— 
’ ’ 73143 


0.5112. 


Giving a limiting value in the context of the IFA approximation 

(10) ^ 16,308,004,933,962,338,624 ^ 
31,975,093,599,832,520,109 

Quite curiously, the renormalization procedure gives also a finite limit when m ^ 0, that is 


-D\^) _ _ n ^^7 TD\^) _ _ n AAA td{^) _ _ 

^ 0.0 ~ 3 — g ~ 27 


0 . 296 , 4 ") = ^ 


0.198. 


Thus, we can conclude that = (|)\ as it is expected from the Cantor structure of the 

counter membrane. Table I presents the values of predicted by the IFA for the first eight 
fractal iterations. It appears that take intermediate values of the order of 0.5 to 0.9. A 
monotonic decay is observed for successive prefractal generations as the distance Ly = ml + 1/2 
remains fixed. It is clear that when Ly becomes very large {Ly oc) Bm tends to unity. As 
an empirical value coming now from our IFA approximation when Ly — 0.1 the coefficient B for 
the limiting fractal object approaches the value 0.51 (see Tables I,II). On the other hand, for a 
given fractal iteration as Ly increases the corresponding values of S’s are increasing. Especially, 
it can be easily shown (through the IFA) that lim Bm =(|)^ + 5 X^(|) • From here, one 

can confirm that for large distances (i.e. m ^ oo or equivalently Ly oo) the limiting value 
of Bm of every fractal iteration is equal to unity, as the spatial distribution of the concentration 
becomes roughly uniform for large distances Ly. Indeed, Table II shows the limiting value of the 
coefficients Bm as the distance Ly is increased. More specifically, we observe a monotonic increase 
for the limiting value from 0.5 to unity as Ly becomes very large (of the order of ten, Ly ^ 10). 
As an illustration of the above discussion. Figure 4 depicts the evolution of the coefficients Bm id) 
(predicted through the IFA) with respect to different fractal iterations for various distances m 
between the electrodes. In each case we observe a decay for increasing fractal iterations and a 
saturation of the corresponding coefficient (after the sixth fractal generation). 


IV. IRREGULAR BOUNDARIES: NUMERICAL APPROACH IN THE NEAR FIELD 


The Laplace equation for geometries associated with the four first generations has been solved 
numerically using an in-house developed finite element code. All four geometries were discretized 
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TABLE I: In the table below, we present the values corresponding to the first few fractal generations of 
the coefficient Bm (g) for different distances from the membrane. We clearly observe a saturation after the 
seventh fractal iteration. It is concluded that the limiting value of for the final fractal object is a 
function of the fractal dimension Df of the object. 

Distance (m) bI^'> 

0.1 0.7222 0.5833 0.5304 0.5150 0.5112 0.5103 0.5101 0.5100 

0.3 0.7917 0.7123 0.6892 0.6834 0.6821 0.6818 0.6817 0.6817 

1.0 0.8889 0.8571 0.8493 0.8475 0.8471 0.8471 0.8470 0.8470 

2.0 0.9333 0.9162 0.9122 0.9113 0.9111 0.9111 0.9111 0.9111 

3.0 0.9524 0.9407 0.9380 0.9374 0.9373 0.9372 0.9372 0.9372 


TABLE II: The limiting value of the coefficient Bm{g) predicted by the IFA is listed for various inter- 
electrode distances (m). 

Distance (m) Limiting value 


0.1 

0.5100 

0.2 

0.6178 

0.3 

0.6817 

0.5 

0.7584 

0.7 

0.8042 

1.0 

0.8470 

2.0 

0.9111 

3.0 

0.9372 

10.0 

0.9794 


using six-noded triangular elements with quadratic shape functions with continuity at interele¬ 
ment boundaries. This resulted in a quadratic interpolation of the potential field and a linear 
variation of its gradient inside each element^^. The global entropy production was then computed 
for various boundary conditions by integrating the local entropy production obtained from the field 
and the associated flux (essentially its gradient) by means of relationship (2). This quadrature was 
implemented using a 25-integration points scheme derived by Laursen and Gellert^^. Dirichlet 
boundary conditions with constant concentrations ci and C 2 are imposed to the models. The 
meshes used for the calculations were produced using the GMSH meshing tool^^. A very fine dis¬ 
cretization of the studied domain was intentionally selected in order to obtain sufficient precision 
on the values of the local entropy production and of the functional for the purpose of validation 
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FIG. 4: Shown are the coefficients B{g) (eq.(16)) versus the fractal iteration i for various inter-electrode 
(m) distances. We observe a saturation after the sixth fractal generation for a range of different distances. 

of the approximate coefficients. An element size on the order of 4.5 10“^ for the points of the 
model was selected (for Lx = 1) close to the fractal surface. The number of elements and nodes 
considered for each generation is reported in Table III (the min and max are related to the minimal 
and maximal distances between the electrodes). 

Using the above method one may estimate the factor B{g) with a very good precision (up to 
three significant digits). The numerical coefficients (obtained from the aforementioned numerical 

method) are presented in Tables IV-V while Table VI shows the comparison of the theoretical (IFA) 

(i) 

and the corresponding numerical values for the coefficients These calculations indicate that 

the proposed approximate method is satisfactory for the first few fractal iterations and gets worse 
for increasing generation. Especially, Tables IV-V illustrate the validation of the phenomenological 
relation with the numerical evaluation via Finite Elements, arguably the central result of this 
paper. From a careful inspection of these Tables it is clear that as Ly increases the coefficients 
increases as well which is at least in qualitative agreement with eq.(16) obtained by the IFA. We also 
observe a decrease of the coefficients A’s and S’s as the fractal generation increases. By inspecting 
Tables IV-V we observe that for Ly ^ 0.1 a limiting value 0.65 is approached from the Finite 
Elements method while the value 0.51 is predicted by the IFA (see Table II). The latter indicates 
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the failure of the IFA to be accurate at high fractal generations. However, a qualitative agreement 
between the IFA and the Finite Elements calculations is evident, which means that for the final 
fractal object the coefficients seem to tend to a limiting value independently of the concentration C 2 . 
Commenting now on Table VI which presents the difference between the theoretical and numerical 
values (obtained with Finite Elements) of the coefficients for C 2 = 0.9, it is observed that as Ly 
tends to zero substantial deviations from the corresponding theoretical value appear as a result of 
the divergence of the Laplacian field. Also, for increasing fractal iterations, we see from the Table 
VI that the relative error becomes important (see also below). 

TABLE III: The number of elements and nodes used for the numerical calculations (Finite Elements) for the 
first four fractal generations. The corresponding values are related to the minimal and maximal distances 
between the two electrodes. 


Generation Number of elements (min- 

-max) Number of nodes (min-max) 

1 

118953-318347 

239450-638544 

2 

174147-385141 

350578-772872 

3 

199691-408749 

402816-821238 

4 

219401-428755 

444236-863250 


TABLE IV: The coefficients and (for the first four fractal generations) obtained from the 

respective numerical values of the variational functional and the phenomenological entropy production 

with concentration C 2 = 0.9. 


Ly 

^{2)ph ^{3)ph ^(4.)ph 

^{l)ph J^{2)ph J^{3)ph j^{A)ph 

0.1 

0.7780 0.6883 0.6620 0.6543 

0.7779 0.6881 0.6622 0.6546 

0.3 

08664 0.8276 0.8170 0.8134 

0.8665 0.8276 0.8172 0.8134 

0.5 

0.9066 0.8833 0.8766 0.8733 

0.9064 0.8827 0.8763 0.8732 

1.0 

0.9500 0.9366 0.9333 0.9300 

0.9491 0.9364 0.9333 0.9301 


Another attempt to improve the above procedure is to insert the phenomenological (numerical) 

values for the coefficients and actually obtain a ’’renormalized” value for the higher order 

( 2 ) 

coefficients e.g. B^^ etc, thus improving considerably the predictions of the theory. Combining 
the empirical result for the first generation with a IFA half penetration argument inside the new 
pores for the second fractal generation, one can find 

^ _ 2 ^ ^ 0.834, (32) 
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TABLE V: The coefficients and (for the first four fractal generations) obtained from the 

respective numerical values of the variational functional and the entropy production with 

concentration C 2 = 0.01. 


Ly 

£{l)ph £{2)ph ^i3)ph ^i4)ph 

^il)ph J^{2)ph J^{3)ph J^{4.)ph 

0.1 

0.7779 0.6882 0.6620 0.6544 

0.7776 0.6874 0.6606 0.6520 

0.3 

0.8662 0.8280 0.8174 0.8139 

0.8661 0.8277 0.8166 0.8128 

0.5 

0.9070 0.8829 0.8758 0.8734 

0.9069 0.8827 0.8755 0.8728 

1.0 

0.9491 0.9361 0.9322 0.9308 

0.9490 0.9360 0.9320 0.9306 


TABLE VI: Percentage difference, with concentration C 2 = 0.9 and various distances between the two 
electrodes for the theoretical and the corresponding numerical values of the coefficients . 


Ly 

B(1) 

5 ( 2 ) 

5 ( 3 ) 

5 ( 4 ) 

0.1 

7.2% 

15.2% 

19.9% 

21.3% 

0.3 

8.6% 

13.9% 

15.7% 

15.9% 

0.5 

8.1% 

11.9% 

12.9% 

13.0% 

1.0 

6.4% 

8.5% 

9.0% 

8.9% 


while 


= 0.786 and = 0.639. 


(33) 


Comparing now the above results we can find that differs from by only 5.8%. 

Similarly differs from B^^^ by only 5.7% and differs from B^P by 8.6%. In the 

same way, we can proceed with the remaining orders. Also, one can try various renormalization 
schemes together with the theoretical arguments. The best results are obtained when one uses for 
the prediction of the coefficient of the A; + 1 fractal iteration only the empirical coefficients 

of the previous k iteration and no theoretical argument at all. For instance, for the case of the 
second fractal iteration for distance m — 0.5 one can find 


(2)fren 

^ 0.5 


0.908 - ^ ) 31 + ^0.949 ~ 0.8746, 

O / o 


(34) 


where the superscript ”fren” stands for ’’fully renormalized”. Here, the difference between 
and ^ 0 ^ 3 ^^ is of the order of 1.2%. Figures 5-12 present the spatial dependence of the entropy 
production for the first four fractal iterations, for Ly — 0.3, for the cases of strong (c 2 = 0.1) 
and mild (c 2 = 0.9) diffusion respectively. We observe that the behaviour of the local entropy 
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production depends strongly from C 2 . In particular, for the case of strong diffusion (see Figures 
5,7,9,11) the appearance of a passive zone is observed irrespectively of the fractal generation. In 
contrast, for the case of mild diffusion the above statement breaks down, being indicated by the 
spreading of the local entropy production near the geometrical irregularities (see Figures 6,8,10,12). 



FIG. 5: Spatial dependence of the entropy production for the diffusion cell of Figure 1. It is shown the case 
of strong diffusion with C 2 = 0.1, while the distance between the electrodes is Ly = 0.3. (Computer units.) 
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FIG. 6: Spatial dependence of the entropy production for the diffusion cell of Figure 1. It is shown the case 
of mild diffusion with C 2 = 0.9, while the distance between the electrodes is Ly = 0.3. (Computer units.) 

Let us briefly recall here the argument which relates the coefficients A and B through the active 
zone concept for the case of mild diffusion, i.e. |(c 2 — ci)/ci| = 0(e), e 1, with ci = 0(1). Here, 
from eqs (7) and (8) one can conclude that P ^ e^ and Q ^ e^ while B(g) ^ lim A^'^^^(ci,C 2 ', g) 

C2^Cl 

up to a factor of Dkjc\. The basic idea of this argument is that the profiles of the concentrations 
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FIG. 7: Spatial dependence of the entropy production for the second fractal iteration. It is shown the case 
of strong diffusion with C 2 = 0.1, while the distance between the electrodes is Ly = 0.3. (Computer units.) 


are essentially linear in space in the whole active volume and the nonlinearities are smoothly 
approximated. 

In the opposite situation, we have the strong diffusion. One particular case of strong diffusion, 
which we study here in detail, is the case of the harmonic measure (ci = 1, C 2 ^ 0). Exactly in 
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FIG. 8: Spatial dependence of the entropy production for the second fractal iteration. It is shown the case 
of mild diffusion with C 2 = 0.9, while the distance between the electrodes is = 0.3. (Computer units.) 

this case we obtain a visual interpretation of Makarov’s theorem. The measure is spreading on the 
surface of the working membrane almost uniformly in a length which roughly coincides with the 
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FIG. 9: Spatial dependence of the entropy production for the third fractal generation. It is shown the case 
of strong diffusion with C 2 = 0.1, while the distance between the electrodes is Ly = 0.3. (Computer units.) 
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FIG. 10: Spatial dependence of the entropy production for the third fractal generation. It is shown the case 
of mild diffusion with C 2 = 0.9, while the distance between the electrodes is = 0.3. (Computer units.) 

diameter of the active membrane / + ^ + ^ + / = 3/. In this case, Makarov’s theorem guarantees 
that at least for C 2 > 0.01, A = B holds (see also Tables IV-V). We thus obtain a very important 
result of our real space renormalization scheme: Estimating the values of B(g) with the active zone 
approximation (Mandelbrot’s argument) we obtain also the value of A in the whole range of C 2 , 
that is in all cases except from the case where C 2 < 0.01. 
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FIG. 11: Spatial dependence of the entropy production for the fourth fractal generation. It is shown the 
case of strong diffusion with C 2 = 0.1, while the distance between the electrodes is Ly = 0.3. (Computer 
units.) 
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FIG. 12: Spatial dependence of the entropy production for the fourth fractal generation. It is shown the 
case of mild diffusion with C 2 = 0.9 and distance between the electrodes = 0.3. (Computer units.) 

The above remarks indicate that the aforementioned IFA scheme is valid (in the Far and Near 
domain) at least for the first few fractal iterations for both strong and mild diffusion. Indeed, it 
is remarkable that such a simple renormalization scheme gives accurate predictions for the first 
fractal generations and leads at least to qualitative results for higher iterations (here up to the 
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fourth). From our numerical calculations it turns out that the quality of the IFA prediction gets 
worse as the fractal generation increases and/or as the inter-electrode distance tends to zero (see 
Table VI). The basic reason behind this result is that the IFA is a coarse-graining approximation 
which neglects the fine microscopic details of the fractal objects. This is exactly the way to the 
passive zone approximation. Thus, it is a challenging open question whether the above conclusions 
can be straightforwardly generalized for even higher fractal iterations, where the IFA approxima¬ 
tion is expected to fail. This constitutes a heavy computational task especially for high enough 
fractals iterations. In particular, three questions remain to be checked: a) to what extent the fully 
renormalized scheme is valid for higher fractal iterations, b) whether the limiting value predicted 
by the IFA approximation holds as a qualitative argument and c) the investigation of a better 
renormalization scheme that takes into account the respective correlations in the distribution of 
the harmonic measure especially for m ^ 0. 

V. CONCLUSIONS AND DISCUSSION 

The transport associated to molecular transfer of neutral particles from a diffusional field to an 
irregular membrane is examined from the standpoint of irreversible thermodynamics. Combining 
non-equilibrium thermodynamics and an argument developed initially by Mandelbrot et al.^^, 
and then by Makarov et al.^^’^^, one is led to probe and predict the entropy production and the 
variational functional of the transfer. Detailed predictions are given for the whole range of distances 
from the Far field to the Near field corroborated by finite differences^’^^’^^ and Finite Elements 
calculations. 

The central idea of the active zone concept^ is that on the grounds of the Laplacian transfer 
on a membrane, certain parts of the membrane remain passive and all thermodynamic properties 
of the membrane can be estimated with a very good approximation by the contribution of the 
’’active parts”. In the context of the active zone introduced by spectroscopic impedance, a system 
with irregular boundaries possesses, in the Makarov regime where the reactive system is ’’almost 
entirely diffusive”, the same impedance as a smaller system with regular boundary conditions. This 
Lapp.Mak ’’apparent, spectroscopic length”, is not the same as the quantity Lapp^entr — B{g)Lx but 
the conclusion remains the same: an entirely diffusive system with irregular boundary conditions 
possesses the same entropy production as a system with smaller diameter and regular (flat) bound¬ 
ary conditions. This can be seen as an interpretation of Makarov’s theorem in terms of irreversible 
thermodynamics. Do these remarks apply to the ’’Near field” ? In the Near field the original Man- 


18 



delbrot’s argument boils down and the notion of spectroscopic impedance is relaxed. However, as 
it is shown here one does not expect ’’new physics” and irreversible thermodynamics is still valid. 
Our remarks, are corroborated from the evaluation of the local entropy production around the von 
Koch curve in the Near field. 

To be concrete, the entropy production and the variational functional of Laplacian transfer for 
geometries associated with the first few fractal iterations leading eventually to a fractal generator 
(von Koch curve) are studied. We first recover and recapitulate results for the Far field. Then, 
we turn to the case of the ’’Near field”, a situation not studied before in the literature, for both 
strong and mild diffusion. Detailed charts of the local entropy production around the first four 
iterations are given, and the global entropy production is investigated carefully. From these studies 
it is clear that the Laplacian fields are in a sense a ’’smooth measure” around irregular objects. 
Moreover, the analytic procedure which allows one to approximate the entropy production seems 
to be quite powerful and straightforward. Although the spatial distribution of the local entropy 
production changes greatly from the mild to strong diffusion, the IFA analysis and the global 
entropy production obey the same rules in the Far and in the Near field. A fair evaluation of 
the IFA arguments implies that the difference between the predicted and the numerically obtained 
values increases with the fractal iterations. An interpretation of the above statement is that the 
IFA is a coarse-graining approximation which neglects the fine microscopic details of the fractal 
objects (passive zone approximation). More specifically, in the framework of our IFA renormalized 
S’s for the first three generations differ by 5 — 15% from the numerical reference values, and the 
fully renormalized B’s coincide up to 2% difference to the numerical values. As already pointed 
out, the generalization of these results for higher fractal iterations constitutes an interesting open 
problem. Another question that arises from our investigation is whether the predicted by the IFA 
limiting value of the coefficient B for high fractal generations holds at least qualitatively. This 
task requires heavy computations for the Laplacian transfer of high fractal objects and has to be 
performed in the context of a supercomputer. Finally, as we have also mentioned in the past, our 
predictions on such dissipation trends should be amenable to experimental testing in a suitably 
constructed diffusion cell. 
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